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ABSTRACT 


Experimental  tables  of  availabilities  at  time  t are 
obtained  for  a device  whose  performance  is  described  by  an 
alternating  renewal  process  with  a finite  number  of  failure- 
renewal  cycles,  corresponding  to  having  a finite  spares 
backlog.  Failure  and  repair  rates  are  assumed  to  be  constant, 
and  attention  is  restricted  to  cases  in  which  the  repair  rate 
is  larger  than  the  failure  rate. 
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I.  INTRODUCTION 


The  most  commonly  encountered  working  definition  of  the 
"availability"  of  a device, 

MTTF 

(1.1)  Availability  * 

MTTF  + MTTR 

represents  the  long-term  or  steady-state  probability  that 
the  device  will  be  found  in  an  "up"  or  functioning  condition 
when  two  specific  conditions  are  satisfied.  One  condition  is 
that  there  is  an  alternation  of  failure  and  repair  cycles 
in  which  times  to  failure  and  times  to  repair  are  independent 
realizations  from  some  failure  and  repair  distributions 
satisfying  minimal  regularity  conditions.  The  second,  and 
here  most  important  condition,  is  that  the  alternation  of 
failure  and  repair  continues  indefinitely,  so  that  the 
performance  of  the  device  is  described  by  a standard  alternating 
renewal  process. 

For  many  equipments,  the  second  condition  cited  above 
implies  access  to  an  infinite  backlog  of  spares.  In  many 
operational  contexts  this  sort  of  spares  support  cannot  be 
realistically  assumed. 

* If  availability  is  considered  to  derive  from  an  alter- 

nating renewal  process  with  a finite  number  of  cycles, 
corresponding  to  a finite  backlog  of  spares,  then  expressions 
for  availability  become  complex  as  compared  with  equation  (1.1). 
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This  thesis  is  devoted  to  computational  experiments  with 
some  "finite  spares"  availability  expressions.  The  end 
objective  of  such  experiments  is  to  be  able  to  determine  the 
circumstances  in  which  equation  (1.1)  furnishes  an  adequate 
approximation,  or  alternately  to  be  able  to  provide  compu- 
tationally feasible  alternatives  to  its  use. 
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II.  MATHEMATICAL  MODEL 

The  mathematical  model  on  which  the  usual  expressions 
for  the  availability  of  a device  are  based  is  an  alternating 
renewal  process  with  an  infinite  number  of  failure-repair 
cycles.  In  situations  where  repair  requires  replacement  of 
a failured  device  by  a spare,  this  corresponds  to  having  an 
infinite  number  of  spares.  The  model  studied  here  is  modi- 
fied to  allow  only  a finite  number  of  failure-repair  cycles, 
corresponding  to  having  a finite  number  of  spares. 

The  simplest  assumptions  about  failure  and  repair  times 
are  made;  failure  rates  are  constant,  and  repair  rates  are 
constant.  Only  those  processes  that  begin  with  a functioning 
device  installed  aTe  considered. 

In  greater  detail,  the  failure -repair  process  considered 
is  as  shown  in  figure  2.1, 


Figure  2.1  Failure -repair  process. 


where  n is  the  number  of  spares,  Tg  is  the  time  to  failure 

for  the  original  device,  T^,  T2, Tn  are  the  times  to 

failure  for  the  n spares,  and  SQ,  S^, Sn-1  are  t*ie 

times  to  replace  the  original  device  and  the  first  n-1  spares 

It  is  assumed  that  TQ,  SQ,  Tn-1*  Sn-1*  Tn  are  inde" 

pendent  random  variables,  and  that  TQ, Tn  are 

exponentially  distributed  with  failure  rate  X,  while 

SQ, Sn-1  are  exponentially  distributed  with  repair 

rate  n. 

The  availability  at  time  t,  A^Ct) , of  the  original 
device,  supported  by  its  backlog  of  n spares,  is  the 
probability  that  the  process  shown  in  Figure  2.1  is  in  an 
"up"  condition  at  time  t;  i.e.,  that  at  time  t either  the 
original  device  or  one  of  its  spares  is  installed  and  still 
functioning. 

The  increment  in  availability  at  time  t due  to  the 
ktlx  spare,  I^Ct)  is  defined  by 


(2.1) 


Ik(t)  - A.  (t)  - A.  ,(t)  k - 1, 


Before  proceeding  to  a derivation  of  expressions  for 
In(t)  and  An(t)  in  a general  case,  two  special  cases  are 
considered;  repair  rate  n equal  to  infinity,  and  repair 
rate  n equal  to  failure  rate  X.  These  are  boundary  cases 
for  the  cases  of  likely  practical  interest,  in  which  it  is 
reasonable  to  expect  that  repair  rate  will  exceed  failure 
rate. 


In  any  case,  AQ(t) , availability  at  time  t with  no 
spares  is  given  by 


(2.2) 


At 


A„(t)  - P[T  > t]  - e t > 0. 


In  the  following  sections,  it  will  be  convenient  to  let 


(2.3) 


Un  ■ T0  * 


* T , 
n* 


vn  ‘ So  * 


* V 


W • U . V 
n n n 


(Sn  + Tn)  * 


+ (S  + T ) 
v n n' 


A.  REPAIR  RATE  EQUAL  TO  INFINITY 

The  simplest  case  is  the  one  in  which  no  time  is  required 
to  repair  a failed  unit,  provided  a spare  unit  is  available. 
In  this  case  the  contribution  of  the  first  spare  is 


(2.4) 


:.(t)  «f  P[TX  > t-s  | UQ  - s]  fjj  (s)ds 

Jo  0 


where 


fri  (s)  * Ae"Xs,  s > 0 , 
u0 


is  the  gamma  { 1 , A > density.  Thus 
(2.5)  I,(t)  - | e”XCt's)Ae'Xs  ds 


l(t)  - r < 

Jo 
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I I 


, -Xs 
Xe  t. 


and  the  availability  of  a system  having  one  spare  is 


(2.6) 


Ax(t)  - A0(t)  + Ix(t) 


= e"Xt(l  + Xt) 


The  contribution  of  the  second  spare  is 


(2 . 7)  I2(t)  » P[T2  > t-s|U1  = s]fy  (s)  ds  , 


where 


,2  -Xs 

fu.  ■ ~rm  . s i 0 > 


is  the  gamma  { 2 , X } density.  Thus 


(2.8) 


I2(t) 


,-X(t-s)  XZ5e'XS  ds 


(Xt)  2e"Xt  yt  ' 


and  the  availability  of  a system  having  two  spares  is 


(2.9) 


A2(t)  * e”Xt (1  + Xt  + ^“)  • 


Generally,  the  contribution  of  the  n spare  is 


(2.10)  In(t) 


PtTn  > t’slUn-l 


0 


Slfy  (S)  dS  , 
n-1 


where 


,nn-l  -\s 

Un-1(5)  ' 5 - ° ’ 


Is  the  gamma  {n,X}  density.  Thus 


(2.11) 


In(t) 


(t-s)  , n n-1  -Xs 


X s e 

iw 


ds 


n -Xt  1 

e — , , 


and  the  availability  of  a system  having  n spares  is 
(2.12)  An(t)  = e'Xt(l  + Xt  + + + ) . 

Availabilities  An(t)  obtained  from  equation  (2.12)  are 
shown  in  table  1.  In  this  table  n represents  the  number  of 
spares,  and  the  contribution  In(t)  of  the  n spare  can  be 
found  by  subtracting  An_^(t)  from  An(t).  The  availabilities 
shown  in  the  last  column  are  for  an  alternating  renewal  pro- 
cess with  an  infinite  number  of  failure-repair  cycles.  This 
corresponds  to  having  an  infinite  number  of  spares. 

B.  REPAIR  RATE  EQUAL  TO  FAILURE  RATE 

The  failure-repair  process  considered  is  that  in  which 
failure  times  and  repair  times  are  exponentially  distributed 
with  equal  rates,  i.e.,  X = q. 
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Table  1 A (t)  for  the  case  nt 


In  this  case,  the  contribution  of  the  first  spare  is 


(2.13)  Ix(t)  « P(T1  > t-s | WQ  = s]fw  (s)  ds  , 


where 


*2  -As 

fw„Cs)  * Hi)  • 5 - 0 • 


is  the  gamma  f 2 , X > density.  Thus 


(2.14)  I.(t) 


■A(t-s)  \2se'Xs 


,2  -At  + 2 
a e t 

~TJTT~ 


and  the  availability  of  a system  having  one  spare  is 


A,(t)  * A_  (t)  + i (t) 


e'Xt  (1  * . 


Generally,  the  contribution  of  the  nth  spare  is 


(2.15)  I 


where 


.(«  ' f PtTn 

J 0 


> t_slWn-l  " slfw  ds  » 


, 2n  2n-l  - As 

£wn-1(s)  * — ran) ’ s - 0 ’ 
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is  the  gamma  {2n,X}  density.  Thus 


(2.16)  Ia(t) 


-X (t-s)  Xe  fl_,2n-l  ^ 

e TTTnr  C J dS 


A2VAt  r 2n-l  ,s 
-nr--  vv  i s as 


'Xt)2ne"Xt 


and  the  availability  of  a system  having  n spares  is 


2 

(2.17)  An(t)  - e'Xt  (1  * £££!_  + 


* (U)2" 
(2nJ ! ) 


C.  REPAIR  RATE  GREATER  THAN  FAILURE  RATE 

The  failure-repair  process  considered  in  this  section 
is  visualized  as  one  in  which  the  repair  rate  is  greater 
than  the  failure  rate.  This  influences  the  format  in  which 
the  results  are  displayed. 

In  this  case,  the  contribution  of  the  nth  spare  is 


(2.18)  In(t) 


f ptT„ 

J 0 


> t-s I W_ 


s]  fw  (s)  ds  , 


where 


ul if  (s)  * / fy  (s)  fjr  (t-s)  ds  , 

n-1  J0  vn-l  n-1 


Table  2 A (t)  for  the  case  nt 


and 


fv„  <s>  ■ rW  sn'VnS  ■ 5 i 0 • 


is  the  gamma  {n,n}  density. 


while 


V , Ct‘,J  ‘ TTn T 


(t-s)n_1  e'X^t's-’  , t-s  > 0 , 


is  the  gamma  {n,X>  density.  Thus 


(2.19)  In(t) 


I mr  “"-Vnu  (.-u)1 


rw 


. e ^t_S^  du  ds  • 


Inverting  the  order  of  integration,  equation  (2.19) 


becomes 


t r t 


(2.20)  In(t) 


0 '■'u 


-x(t-s)  nn  ,.n-l  -nu  xn  ^ 

rw  e rinT  (s 


s-x(s-u)  ds 


-xt  nnxn 


(r(n) 


n rt  rt 

urx-le-(n-X)u  ( (s.u)n”l 

0 Ju 


Let  v * s - u.  Then 


s 


Thus  equation  (2.20)  reduces  to 


(2.21)  In(t)  - e 


-xt  nn*n 


. f uP-^t-u)11 

J 0 


e-(n-x)u  du 


and  the  availability  of  a system  having  n spares  is 


(2.22)  An(t)  - A0(t)  + Ix(t)  + 


* • 


Note  than  when  n » X,  equation  (2.21)  reduces  to 


equation  (2.16),  since  then 


XnCt)  * e 


-Xt  X‘ 


__  uJ1“1(t-u/)n  du  . 


Let  u .»  tv.  Then 


r t 

J 0 


(tv)n’1{t(l-v)}11  tdv 


■fVw 

Jo 


c2n  w 
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III.  APPROXIMATION  OF  THE  MATHEMATICAL  MODEL 


In  section  II,  we  derived  a mathematical  model  for 
availability. 

In  this  section,  we  discuss  methods  of  approximation  to 
obtain  numerical  values  of  availability. 

A.  EXPONENTIAL  EXPANSION  APPROXIMATION 

The  integral  in  equation  (2.21)  can  be  approximated  by 
expanding  its  exponential  term,  i.e., 

(3.1)  e'Cn-X)  , x _ (n_X)u  + (n-.j)2  ul  . 

Thus  the  integral  becomes 


du 


t 

un(t-a)n  du 


(n-x) 

TT~ 


r* 


...n+l,*  <>11  , 

u (t-u.)  du 


} du 


We  know  that, 


19 


t 


I 


1 ..n-l^  .^n  J..  _ „2n  r(n)r(n+l) 

J u Ct_u)  du  t T'C2n+T) 


Thus  equation  (3.2)  becomes 


(3.3) 


^ -n.2n+l  r(n+l)r(n+l) 
Cn_;gt  rczn+T) — 


* (n-x)2t2n+2  r (n+2) r (n+i) 
21 HTn+TJ 


e"Xt(nt)n(Xt)n  r T(n)  ^ T(n+1) 

— mn {r(2nTrj  • Cnt  xt)  wit 


. (nt-xt) L r (n+2) 

21 r (2n+TF 


Computational  experiments,  with  the  approximation 
represented  by  equation  (3.3)  have  indicated  unsatisfactory 
convergence  behavior  when  (nt-Xt)  is  large,  a case  of  some 
practical  interest,  and  so  this  approach  was  not  pursued. 

B.  SIMPSON'S  RUL2  APPROXIMATION 

The  integral  in  equation  (2.21)  can  be  approximated  by 
using  Simpson's  rule. 


I 
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Let  v 


Then 


u 

t * 


(3.4) 


f u~  WC--H.  du 

J 0 


,n-l 


(tv)  .{t(l-v)}n  e-(n"X)tVt  dv 


r i 


.2n 


rn-l(1.v)n  e- (nt-xt)v  dy  . 


Thus  equation  (2.21)  becomes 


(3.5)  In(t)  - e'xt  5 iSj  1' Jnvir 


vn(l-v)n  e-Cnt-Xt),r  dv 


Now  Simpson's  rule  can  be  applied  to  the  integral  in 
equation  (3.5)  to  obtain  numerical  values  of  availability. 
Simpson's  rule  as  applied  is 

C3.6)  f f(x)dx  - j(y1*4y2'-2y3*4y4*....*2ym.l*«ym*ym<.1) 


a 


where 


b-a 

m 


and 


Xx  - a,  X2  - a+h,  Xm+1  - a+mh  - b, 
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IV.  TABLE 


In  this  section  experimental  tables  of  availabilities 

A ft)  are  shown,  which  were  obtained  from  the  mathematical 
n 

model  evaluated  using  Simpson's  rule,  i.e.,  using  equation  (3.6) 
to  evaluate  equation  (3.5)  with  m * 500,  ■ 0.0001  and 

X501  *0.9999.  The  reason  for  choosing  m - 500  is  that  compu- 
tational experiments  with  choices  of  m greater  than  or  equal 
to  500  gave  a stable  and  accurate  result. 

In  these  tables  n represents  the  number  of  spares, 
and  the  contribution  of  the  nth  spare  In(t)  can  be  found  by 
subtracting  An_^Ct)  from  An(t) . 

The  availabilities  shown  in  this  section.  Table  3-12, 
are  for  cases  in  which  nt  > Xt. 

The  availabilities  shown  in  the  last  column  are  for  an 
alternating  renewal  process  with  an  infinite  number  of  failure- 
repair  cycles . This  corresponds  to  having  an  infinite  number 
of  spares. 
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Table  3 A (t)  for  nt 


Table  5 A ft)  for  nt 


Table  6 A (t)  for  nt 


Table  8 A (t)  for  nt  a 120 


Table  9 Aft)  for  nt  = 140 
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Table  11  A (t)  for  nt  = 180 
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Table  12  A (t)  for  nt  * 200 


r.  SUMMARY  AND  CONCLUSIONS 


l - 

Certain  computational  approaches  were  tried  for  obtain- 
ing availabilities  for  a device  supported  by  only  a finite 
backlog  of  spares,  using  the  simple  assumptions  that  failure 
and  repair  rates  are  constant.  Real  failure  and  repair  distri 
butions  may  be  more  complex,  but  the  case  considered  is  a 
good  case  for  initial  computational  experiments. 

Of  the  two  approaches  tried,  neither  proved  entirely 
satisfactory  in  obtaining  availabilities  in  a way  that  is 
fast  and  suitable  for  use  with  small-scale  computational 
facilities,  e.g..  hand-held  calculators.  Also  neither  was 
effective,  over  the  entire  range  of  failure  rate  and  repair 
rate  combinations  that  might  be  of  interest. 

Since  an  easily  used,  readily  accessible, way  to  assess 
the  impact  of  finite  spares  backlogs  on  availability  is 
desirable  in  many  mission  planning  contexts,  further  computa- 
tional approaches  should  be  tried. 

The  tables  presented  in  section  IV  give  availability 
values  with  which  the  result  of  such  experiments  can  be 
compared. 
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